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ABSTRACT 

Both Lord (1967) and Cronbach and Furby (1970) have shown 
that there is no knowably correct method of analysis- for data 
from a nonequivalent control group design unless some assump- 
tions are made. The assumption made in this paper is that of 
a particular class of continuous growth models . It is shown 

that all -of the traditional methods of data analysis are 
v. 

incorrect under much of this class of growth models. New 
methods of data analysis are then developed based on maximum 
likelihood estimation, jackknifing, and numerical analysis 
techniques. { 
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This paper is concerned with estimation and hypothesis 
testing of treatment effects in nonequivalent control group 
designs (Campbell, 1969; Campbell & Stanley, 1966). A nonequi- 
valent control group design is a quasi-experimental design in 
which the groups are formed by some method other than random 
assignment. Both Lord (1967) and Cronbach and Furby (1970) 
have argued that unless some assumptions are made there is no 
knowably correct method for the estimation and hypothesis 
testing of treatment effects^or these designs. In this paper 
it is assumed that, in the absence of treatment effects, 
natural growth conforms to a- particular class of continuous 
growth models. 

This class of continuous growth models can be expressed 
symbolically as, for all times t, ^ 

Y*j(t) = g. (t)-Y* j (t 1 ) + h. (t) + a. (t) 

and ' (1) 

^(t) = Y* (t) + e.j(t) ; j = l,2,...,J 

where J Is the number of groups in the design; 

g^ (t) and h_. (t) are continuous functions; 

a., (t) represents the population treatment effect 

for the jth group; 

t is an arbitrary time point; 
j 

ri n d ^ 

Y^j (t) , Y^j (t) , and e^j (t) represent the .true scores, 

observed scores , and errors of measurement, respectively , 
for the iift individual in the jth group, on the measure 
of interest. 

/ 



2 

Further assumptions are: ^ 
(1) Classical measurement theory holds • That is, for each 
time t, Yj(t) and e.. (t) are uncorrelated and Efe^tt)) = 0 • 

and (2) Treatment effects are additive. 
Notice that for this class of growth models, the correlation 
within each group between true scores at any two points in time 
is +1 . Further, notice that any number of groups are allowed 
under this class of growth models and that g^ (t) • Y^.. (t) + h^ (t) 

represents natural growth for these models. Finally, notice 
that thl treatment effects, as defined by the a_.(t) f s in the 

system of equations (1) are not the same as the usual defini- 
tion of treatment effects. Let a(t) be the grand mean of the 
<Xj(t) v s. The usual definition of a treatment effect is given 

by ctj (t) - ct(t). Although it might be argued that this clarss 

of growth models is extremely restrictive, the authors (Blumberg 
& Porter, 1981) have shown that, in fact, a wide variety of 
natural growth patterns are included in this class (e.g. 
differential linear growth, exponential growth, and logistic 
growth) . 

The appropriateness, under the class of growth models . 
described by the system of equations (1) of currently avail- 
able methods of data analysis has been explored by Blumberg 

(in progress) . She showed that ANOVA of Index of Response 

J 

with K = g(t) (Cox, 1358), ANOVA of Standardized Change Scores- 
with reliability correction (Kenny & Cohen, 1980)' and ANOVA 

ERjC r 6 



of True Residual Gain Scores (Cronbach & Furby, 1970) yield 
unbiased estimates of differences in treatment effects and 
correctly test whether or not these differences &re statis- 
tically significant when the system of equations (1) reduces * 
to * 

Y* (t) = g(t)-Y* (t ) + h(t) + a.(t) 

Y..(t) = Y*.(t) + e ij( t) 

and * * 

0 Y (t) = a y (t) for all j 

j t 

i 

where (t) represents the standard deviation of the true 

j 

scores for the jth group. 

Further, ANOVA of Raw Residual Gain Scores (Cronbach & Furby, 

1970) and ANCOVA provide unbiased .estimates of differences t 

in '.treatment effects and correctly test for nonzero differences 

v 

when ; 

Y*^t) = g(t) -Y*.. (t L ) + h(t) + a.(t) 
0* (t) = 0*{t) for ail j 



apd 



no errors of measurement are present in the data. 
Other special situations were also found where one or 
more of the existing analysis strategies provide unbiased 
estimates of treatment effects and appropriate significance 



testing procedures. Nevertheless, none of the conventional 
.approaches to analyzing data from nonequivalent control group 
designs -are appropriate when natural growth .conforms to the 
system of equations (1) unless the gj(t) f s are all equal and 

the hj(t) f s are all equal. It should be pointed out here that 

the method of Empirical Bayes Estimation (Bryk, Strenio, & 
Weisberg, 1980; Strenio, Weisberg, & Bryk, in press) seems to 
be a very promising approach for the analysis of data arising 
from the application of nonequivalent control group designs 
'under a wide variety of conditions^. The method of Empirical 
Bayes Estimation, however, requires that the variance-covariance 
matrix of the true scores be non-singular. When the system of 
equations (1)> is assumed, this variance-covariance matrix becomes 
singular and hence the method of Empirical Bayes Estimation can- 
not be used for the class of continuous growth models being 

considered *here . 

Hence," there is a need for new methods of data analysis 

which can provide point and interval estimates of and appro- 
priate hypothesis testing procedure^\for treatment effects 
under the class of continuous growth lnodels considered here. 
It is, however, not possible to develop methods which will 
work for the entire class of continuous growth models since 
the hj (t) and otj (t) terms are confounded. Thus, wften consider- 
ing methods for estimating and testing treatment effects, assuminc 

» 

natural growth as defined in the system of equations (1), it is 
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necessary to distinguish between the following possibilities 
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the hj (t) terms : $ 
(a) the exact natures of * the hj(t)'s are known (i.e., 



h x (t) = 3-t, h 2 (t) = 4-t 2 + 3;t h , ... , hj(t) - log 10 (5-t 3 + 1)); 

(b) the functional forms of the h_.(t)'s are known (i.e., 

h x (t) = kj't, h 2 (t) = k 2 -t 2 + kg-t*, ... , hj(t) = log^fk^t 3 + 1) / 

where k^, k 2 f k^/ and k^ $re unknown real-valued constants); 
and 

(c) for each time tj the h_,(t)'s are equal to som'e common 
value, say h(t) (i.e., for each t, h^(t) = h 2 (t) = ... = hj-(t) 
= h(t)) . 

The discussion of methods of data analysis for each of the above 

possibilities can be further broken down into six cases according * 
♦ 

to whether the exact natures of the g^ttj's are known, the 

functional forms of the g^ttj's are known, or nothing is ki^own 
about the g\ (t) 1 s ^nd according to whether or not errors of 
measurement are present (see Figure 1) . 

The remainder of this p^per will present analysis procedures 
to be used for cases 2, 4, .6, 8, 10, and 12 (see Figure 1) , since 
these are the cases which are likely to occur for data arising 
from educational research settings. Cases 14, 16, and 18 are 
also of interest, but the procedures to be used in these cases 
are sufficiently different that they merit a separate piece. • 
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Figure 1 

; Subcases of the general growth model 
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Procedures have also been developed for the cases where no 
errors of measurement are present. The interested reader is 
referred to^Blumberg (in progress) for the procedures for the 
cases £Ot to be discussed, here . 

V 

Overview of the Procedures 

Each of the procedures to be discussed is a variation of 

a several stage method. At the first stage, estimates are 
t 

obtained for the unknown constants (i.e., parameters) in the 

functional form expressions for the g.(t)'s and/or h.(t)*s. 

3 3 

The method used to obtalin the estimates* of the parameters in 
the functional forms* expressions for the g. (t) 's and/or h. (t) 

is the sarnie whether one is estimating the parameters for the 
g i (t)'3^ for the h. (t)'s, or for both the g.(t)'s and h.(t)'s. . 

\ 

Hence, the discussion" of the first stage will be done in 
general. The second stage concerns the "estimation of treatment 
effects (i.e., the a., (t) 's in equation (1)), given the estimates 

^obtained at the first stage and the particulars of the case of . 
interest. Hence fox; this second stage, each of the cases must 
be discuSseU separately. The third stage concerns methods for • 
interval estimation and hypothesis testing of both treatment 
effects ajnd differences in treatment effects. These methods < 
are the same for" ^11 the eases, once estimates qf the a. (t)'s 

have been obtained. Therefore, the discussion of the third 
stage will be dorie' for all § cases simultaneously. - . 
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Stage 1: Estimation Of the -g^(t) f s -and b^.(t)'s 

i 

The methods developed at stage 1 depend upon knowledge of 
each individual 1 s -scores on several pretests, DenotV-the times 
of the pretests for the jth group "by t- , t- , t- , • • • , t , v 

where p.. is the number of pretests . Hence, for the jth group, 

the system of equations (1) can be rewritten for the pretest times 
(remembering that otj(t)EOfor all pretest times) as 

r 

VV* - v t *. ) - Y ij (t i. ) + h j (t k\> 

3- 3 3 3 

and v J (2) 



Y ij (t k. } " Y ij (t k. } + e ij (t k. } ? k=1 ' 2 Pj ' 



3 3 . 3 

- i=l , 2 , . . . , N . 

3 



( 



where is the number of individuals in group 

Stage 1/is divided into two substages. At the first 

substage, estimates of the g . (tL )'s and h.(t, )'s are ob- 

3 k. j j 

tained. At the second substage, the estimates of the g.(t )'s 

3 

and h.(t, ) f s are used to obtain estimates of the unknown con- 

3 j 

stants in the functional form expressions of the g.(t)'s 

r 3 . 

and/or h . (t) 1 s* 
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For simplicity, the j subscript will be 'dropped for the 
remainder of the discussion of stage 1, since the estimation of 
the unknown constants in the functional form expressions for 
the g.(t) f s an,d h.(t)'s is done separately fo^ each of the J 

groups. Dropping the j subscript from equation (2), ' 

Y*(t k ) * g(t k ):Y^(t 1 ) + h(t k ) (3) 



and 



Y i ( V - W + W ' k=1 < 2 



(4) 



Equation (3) represents .a lihear structural (Madansky, 1959; 
Moran, 1971) or functional relation (DeGracie & Fuller, 1972; 

Lindley, 1947), where ^(t^) and h(t^) are the slope and 

* * * * 

(t^) -intercept , respectively, of the Y^(t^) on Y ^( t i ) 

regression line. There are several approaches that have been 
developed for finding estimates of the g(t k )'s and hft^h's. 

The ^£proach that is to be us^d here is maximum likelihood 
estimation* 

m 4 

Estimation of the g^(t^) 1 s and hCt^-'s 

For a design with p pretests the system of equations 
represented by equations (3) and (4) can be described pictorially 
as in Figure 2. ' 
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!> . Y*(t 2 )^ ^ e ^ 2 ) 




Figure 2 

* » 

Pictorial representation of the structural relation. 
The system can be^ written in vector form as 

and 

" W ■ VV + e i ( V . • • 



' . . (5) 



(6) 
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Maximum likelihood requires expressions for the means, variances 
and covariances of the observed variables under considera- 
tion. Taking the mean on both sides of equation (5) gives 

- U Y (t k ) = ^V'^V + h( V (7) 




where U y (t k ) is the v^?tor of means. The variance of Y(t^) 



is given by 



ajttj) = [a*(t L ) ] 2 + a^(tj) • (8) 

2 

where (t) represents the variance of the errors of measure- 
ment at time t. The variance of ^(t^) is given by 

ayCt k ) = [g(t k )] 2 -[a*(t x )] 2 + a 2 (t k ) (9) 



for k=l,2, ... , p . The covarinnce of ^(t^) an ^ Y(t^) ,is 
.given by 



a Y(-t k )Y.( tl ) -'« ( V I0 i (t i! 1? • " (10) 

The covariance of Y(t k ) and Y(t k ,), where k,k'.= 2,3, ... , p 



and k f k 1 , is given by 



4 



°Y(tjy(t,) =.g(t k )-g(t k ,)-[a; ( t 1 )] 2 . 
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• Cdnsider the system of equations (7) ; (8) , (9) (10) , and 
(11) as one large system of equations* The maximum likeli- 
hood approach requires that this system be identifiable, A 
system of equations is said to be identifiable if and only 
if each of the unknown paramethers 'on the right hand side 

i 

of the system can be expressed in terms -of the unknown 
parameters on the left hand side of the system. For the 
models of continuous growth being considered here, the 
unknown parameters on the left hand side are y^ft^, Pyftj)/ 

^y(tp) an ^ the P - parameters in the variance- 

covariance matrix of Y(t n ), Y(t n ), • , Y(t )♦ The unknown 

12 i p 

i 



parameters on the riglft hand side are g(t2)/ g(t^)/ 



g(t >, h<t 2 ), h(t 3 ) ; ... , h(t ) r [o*^)] 2 , o*(t x ), ^(t 2 ), 

2 2 
. .* , a (t and a (t ) . The system of equations (7) to 

; e p-»r e p 

(11) is identifiable if and only if p > 3. When p >_ 3, the 

LISREL program (Joreskog & Sorbom, 1978) can be used to 

obtain the maximum likfelihood estimates of the parameters qh 

the right hand side. When p = 2, additional assumptions must 

be made in order that*the system of equations is identifiable 

\ 

Estimation of the Unknown Constants 

» 

Once the'' estimates of the g(t^) ! s and h(t^) ! s are ob- 

■ ; ' ^ — 

tained, estiflatesyof the unknown constants fin the general 



\ 
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functional form expressions for g(t) artd h(t) can then be 
computed. For simplicity, the discussion here will be in 
terms of g(t). The method to be used for h(t) is analogous. 

^ The method for determining estimates for the unknown 
constants in the expression for g(t) is dependent upon the 
functional form of g(t). For sake of illustration, consider 
the following two examples, 

. 3 

a polynomial form: g(t) =1+ E c • (t - t ) (12) 

d=i d 1 

ana i . :. 

an exponential form: g(t) = b # c (t ~ t l ) + (1 - b) . (13) 

Polynomial Form 4 

, The object of the method to be described is to determine 
the values of c^, c^ r and c^ in the polynomial from knowledge * 

of the gTt^Vs, where the g7t^)'s are the estimates of the 

g(t^)'s obtained using the maximum likelihood approach. For 

simplicity, assume that p-5, t^ = 0, = 1/ ^ = 4, t^ = 5, and 

t 5 =7. The method to be described will work for any p > 4 

and for any set of values for bhe t^' 5 * Substituting t^ = 0, 

t 2 = 1, t^ = 4, t^ = 5, and,t' 5 = 7 into equation (12) 

» 

yields 
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(14) 
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g(l) = 1 + c 1 + c 2 + c 3 

g(4) = 1 + 4c 1 + x 16c 2 + 64c 3 

g(5) = 1 + 5c 1 + 25c 2 + 125c 3 
and 

g(7) = 1 4 7c. +,49c„ + 343c. . 

\ 

But, the values of g(l), g(4), g(5), and g(7) are unknown. 
Therefore, the* maximum likelihood estimates are used instead, 

Replacing g(l), g(4), g(5), and g(7) by gTT) , g?4) , g?5) , 

and g?7V in the system of equations (14) yields 

/ 



g?l) = 1 + c 1 



and 



1 + c, + c 2 + c^ 



g(4) = 1 + Ac 1 + 16c 2 + 64c 3 



g(5) = 1 + 5c 1 + 25c 2 + 125c 3 



g(7) = 1 + 7c x + 49c 2 + 343c 3 



(15) 



The method of least squares is then used to obtain estimates 
Of c^, and Cg from the system of equations (15) 



• V 



Exponential Form 



The object of the method to be described in this sub- 

section is to determine the values of b atid c in equation (13) 

from .knowledge of the gTtTj's. The method to be described 

K * 
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will work for any p > 3 and for any set of values for the 

t^'s. Recall that the g(t^)'s were already determined at 

the first subrtage. Next, the values of the t^'s and 'gTt^) 1 s 

are substituted into equation (13)*- to yield a new system of 
equations analogous to system (15) . For sake of illustration, 
let p = 4 , t^=0, - If t^ - 4, and - 5. Substituting 

these values into equation (13) yields 
gU) = b-c + (1 - b) 



and 



(4) = b-c 4 + (1 - v b) (16) 



g(5) = b-c 5 + (1 - b) 



Estimates of b and c are then derived from the system of 
equations (16) using a combination of the numerical analysis 
techniques of least squares and the Newton-Raphson method. 

V 

\ 

4 

Stage 2: Point Estimation of Treatment Effects 

The discussibn of the determination of the. point 
estimators of treatment effects will be discussed separately 
for each cask, since the process is slightly different in 
each case. 
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Case 2 estimates 



For case 2 the exact natures of the g.(t)'s and h.(t)'s 

J 3 

are assumed known* Unbiased point estimators of the treat- 

0 



ment effects are given directly by 




a.(t) = Y.(t) - [g (t) -pr . (t x ) + h.(t)] 
Notice that for this case, stage 1 is skipped completely, 
Case 4 estimates 



Case 4 assumes the exact natures of the h.(t)'s are 

3 

known but that only^jJ^^nctional fotms of the gj(t)'s are 



\ 



knovtn. Recall that at stage 1, estimates were obtained for 

the unknown constants in the functional form expressions for 

the g.(t)'s* For) the jth group, let g.(t) denote the 
3 > - 3 

function formed by substituting these estimates of the 
constants into the general functional form expression for 
g.. (t) . Point estimators of treatment effects are then given 

by \ 



'a. (t) = YTTET ' - [§. (t) -Y. (t x ) + h. (t) ] . 

Case 6 estimates 

_ Case 6 assumes that the exact natures of the gj(t)'s 
are known but that only the functional forms of the h_.(t)'s 
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are known. Analogous to case 4, for the jth qrdxxp, form 

7 

a new function, called fL (t) / by ■ substituting the estimates 

of the unknown constants that were obtained in stage' 1 into 
the functional form expression for h.. (t) . Point estimators 

of treatment effects are then given by 



cx^t) = y. (t) - [gjOO-Yjd^ ) + fi^t)] 



Case 8 estimates 



Case 8 assumes that only the functional forms of both 
the g. (t) *s and h. ( are known. For this ".case, point 

estimators of treatment effects are given by 

vAX) = YTTtT - (t) • YTTt^T + fi. (t) ] 

where the §.(t) ! s and fi.(t.) f s are as defined previously. 



\ 



Case 10 estimates 

1 

Case I'O assumes that the^xact natures of the h.(t) ! s 

* 1 3 

>are known and that nothing is kiiown about the g_j(t)'s. j 

Define a new variable, W i ^ (t) - h.. (t) . Then, the 

system of equations (1) can be rewritten as 



22 
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W*..(t) = g.. (t) .W*^ (t 1 ) +a j (t) 



and (17) 

V 

W. . (t) = W. . (t) + e. . (t) . 
ID ID « ID 



Notice that the system of equations (17) is, for each j , a 

* 

linear structural relation with a slope of g.. (t) and a W..(t)- 



intercept of a_.<t).' Hence, for any particular time t, esti- 
mates of g.(t) and a.(t) can be obtained directly by using 
I 3 3 

maximum likelihood techniques, with W.(t, ) as the inde- 

3 *j 

pendent variable and W.. (t) as the dependent variable. 

should be recalled here that in order to compute the maximum 
likelihood estimates in situations where there are only one 
t independent and one dependent variable, it is necessary to 

make additional assumptions. The maximum likelihood estimates 
of the aj(t)'s are considered to be the stage 2 estimates and 



are labelled ct_.(t) 's. 



Case 12 estimates 



>t^ma 



Case 12 assdmes that only the functional forms of the 
h_.(t)'s are known and that nothing is known about the g_.(t)'s. 

In this case a new variable is formed by defining U^(t) = 

Y ij^ ^j(t)/ where fij(t) is defined as in case 6. The 
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method described for case 10 is then used to obtain ~the 
<*.(t)'s by replacing the W..(t)'s o£ case 10 with the 
^.(tj's. 

In developing the point estimators of the treatment 
effects it was assumed that for any particular nonequivalent 
control group design, the known information about the growth 
curves for each of the J groups belonged to the same case. 
Since the determination of the point estimators is done 
separately for e^ch group, insisting that all of the groups 
belong to the same case is overly restrictive. Hence, it . 
will be assumed for the remainder oi this paper that for each 
of the J groups, the known information about the growth curves 
allows the data analyst to place each of the groups into one, 
of the cases 2,4,6,8,10, or 12. Once this is done, the ^.(t) 1 

are derived separately for each group by using the methods 
given in this stage 2 section for the case to Which the ♦ 
group's growth curves belong. 

Bias of the Point Estimators 

The procedures used to determine the -point estimators 
of the treatment effects, involved the use of maximum likeli- 
hood estimation at either stage 1 and/or stage 2, except for 
case 2 type growth curves. It is well known that maximum 
likelihood techniques often lead to biased estimators.: These 
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estimators are, however, usually consistent (Patel, Kapadia, 
& Owen, 1976) and asymptotically efficient (Bickel,& Doksumj, 
1977) . Further, functions Qf consistent estimates are qilso 
usually consistent estimators. The asymptotic properties 
of the point estimators for treatment effects derived here 
for the various cases have not been studied in much detaii^^ 
Further- study of the asymptotic properties is needed. 

Stage 3; Interval Estimation and Hypothesis Testing 

M 

In order to develop interval estimation and hypothesis 

testing procedures it is necessary to have estimates of the 

variances of the a .(t)'s. The method the authors have chosen 
/ 3 

to use to obtain estimates of the desired variances is jackknifing 

1 V 

(Quenouille, 1956) . The method of jackknifing was chosen for 
two reasons. Fir£fc, it is the only practical method of ob- 
taining variances for estimators that are formed by the use 
of a several-stage process such as| those developed here. 
Second, and more important, the use of jackknifing usually 
leads to new point estimators ^hich have a reduction in bias 
over the original point estimators, when the original esti- 
mators are biased. 

k 

The technique of jackknifing begins by drawing a random 
sample from a specified populati^h. Let N denote the number 
of subjects in the sample. The N subjects are then divided 
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. . . . . N ' 

into m disjoint subsets, each of size . Let y'be the, , 

m , 

parameter of interest and ^ be an estimator of y. Further, 

T • * 

let ? be the value of ? when all N subjects are us§d, and 

let f ^ be the value of ? when the subsairiple of size N — N 



m 

where the i th subset has been deleted, is used. 
Next define J^(?) by 

; J^(?) = m-? T - (m-l)? U) . ; £=1,2,..., m 

I 

and' define J {fy by 



Aii estimate of the variance of J(?) is given by (Tukey, 1956) 

m 2 

S J ~ 



m - 1 



Think of m as being fixed. Then as N approaches infinity, 

NT * * 

-jjj- also approaches infinity. Gray and Schucany (1972) have 



shown 



that ^ tJ( ^ ^— is asymptotically distributed 



(as — ► °°) as a random variable with m - 1 degrees of 

m 

freedom. 
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The estimators of interest here are the ct.(t) f s from 

stage 2. To apply the -jackknif ing technique, \frirst divide 
the Nj subjects from the jth group into distinct subsets. 

Next/ define cc^ (t) to be the value of the -estimator, "cT^TtT^ 
yhen the &th subset is deleted from the -sample ,for the j th 

group. O^tce the ft (t) f s have been computed, new estimators 
of the aj(t) ! s are given by ^ 

J(aT(t)) = -p- Z (18) 



/here 



J £ (cc.(t)) = nya.(t) -^f- 1) 'a^ (t) 

Y A 
Interval estimation and hypothesis testing procedures are / 

available by observing that 



/ny [J(o.(t) ) - o. (t) ] 




/ (S J (t)) j ^ 



is asymptotically distributed as £ random variable with a 
Student's t distribution with m.. — 1 degrees of freedom 

; ' I [J 0 (O.(t)) - J(0.(t)) ] 2 . 

where (S^(t)). = * 1 ' = . . 

• J 3 m. - 1 



3 
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An approximate (1- — a) % confidence interval for a_. (t) (and 
hence, an ct-level, test for non2ero a.(t)) is then'given by 



a 1 J- 



j(*.(t)) +. t . (i - ■—) • r— 

i — m.— 1 l / ra . 

J 3 3 



The estimation and hypothesis testing of differences in 

treatment effects are also of interest. A point estimator of 
( 

ct.(t) — ct.,(t), for two groups j and j', can be given by 

J(cu7t)) - J(aj,(t)) where J(oTTTtT) and J(ct_. f (t)) are defined 

by equation (18). An approximate (1 — ct)% test of the hypothe- 
sis of nonzero differences in treatment effects is accomplished 

by performing an ANOVA on the J 0 (ct.(t))'s. Notice that the 

unit of analysis for the 1 ANOVA being performed here is the 
disjoint subsets formed in order to do the jackknifing and 



that the dependent variable is J^(a_.(t)). 



Summary 



/ 
I 



In this paper point estimation, interval estimation, and 
hypothesis testing procedures were developed for both treatment 
effects and differences in treatment "effects when 

Y ij (t) = g j (t) * Y ij (t i . } + h j (t) + a j (t) 



and 



Y. .(t) = Y?.(t) + e. .(t) 
13 13 13 • 
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The only additional assumption needed was that for each of the 
J groups either the exact nature of or the functional form of 
hj(t) was known. Hence, techniques were developed which are 

applicable for a variety of natural growth situations where 
none of the presently available data analysis methods can be 
employed. 
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